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Abstract 

Motivated by Higgs Portal and Hidden Valley models, heavy particle dark matter that 
communicates with the supersymmetric Standard Model via pure Higgs sector interac- 
tions is considered. We show that a thermal relic abundance consistent with the mea- 
sured density of dark matter is possible for masses up to ~ 30 TeV. For dark matter 
masses above ~ 1 TeV, non-perturbative Sommerfeld corrections to the annihilation rate 
are large, and have the potential to greatly affect indirect detection signals. For large 
dark matter masses, the Higgs-dark-matter-sector couplings are large and we show how 
such models may be given a UV completion within the context of so-called "Fat-Higgs" 
models. Higgs Portal dark matter provides an example of an attractive alternative to 
conventional MSSM neutralino dark matter that may evade discovery at the LHC, while 
still being within the reach of current and upcoming indirect detection experiments. 



1 Introduction 



Recently, there has been a surge of interest in models where the Standard Model (SM) or 
the Minimal Supersymmetric Standard Model (MSSM) communicates with a partially 
hidden sector via either Z' or Higgs interactions [H [21 [Sj HI O [6] . These Hidden Valley or 
Higgs Portal models provide a stimulating and consistent alternative to the usual model 
building assumption of a desert above the weak scale. Higgs-sector and Z' interactions 
between the hidden sector and the SM states are special in that they involve gauge- 
invariant operators of dimension do < 4, and thus can be induced by physics at arbitrarily 
high scales with unsuppressed couplings. In the case of a Z' the interactions can either 
occur directly with SM states if they are charged under the f/(l)' or, possibly more 
interestingly, indirectly due to a kinetic- mixing term, eF^F'^^^ between hypercharge 
and the new in which case e, and thus the effective size of the SM-hidden sector 

interaction, can be suppressed [3 [H E]. On the other hand, in the case of the Higgs- 
sector interactions of interest to us here, couplings of the form |ifps^ involving the SM 
or MSSM Higgs states and new SM gauge singlet states can be large, especially in the 
situation where the TeV-scale theory UV-completes not far above the weak scale to a 
strongly interacting theory with light composite states. 

It is interesting to ask whether such models lead to new dark matter candidates with 
qualitatively different phenomenology. In this paper we argue that dark matter commu- 
nicating with a supersymmeterized SM purely via Higgs-sector interactions (the Higgs 
Portal) leads to new and unusual features0 First, as we will show, the thermal relic 
abundance in such a scenario can be consistent with the measured density of dark mat- 
ter for masses as high as ~ 30 TeV, much larger than are usually considered (while also 
being consistent with the upper bound on the mass of thermal relic dark matter derived 
from unitarity [H]). Second, for dark matter masses above ~ 1 TeV non-perturbative 
Sommerfeld corrections [12] to the low-velocity annihilation rate are large. Several au- 
thors have recently recognised the potential importance of these corrections to the dark 
matter relic density calculations |T3l [13 [IS], which lead to enhanced annihilation 
rates in the case of attractive interactions. Even more importantly, as we will argue in 
detail in a companion paper [17] , these corrections have the potential to greatly enhance 
the indirect annihilation signals by factors of up to 10^ beyond those predicted without 
consideration of the Sommerfeld factor, potentially leading to a significant change in the 
optimal search strategy. 

As well as providing examples in which the dark matter particle is beyond the kine- 
matic reach of the Large Hadron Collider (LHC) but is potentially detectable by indirect 
and direct dark matter searches, the models presented here are independently motivated 
by the desire to raise the MSSM upper bound on the lightest Higgs mass, and so relax 
the current tension with the LEP2 Higgs-mass exclusion limit. It is also interesting that 
our models may be given a UV completion in so-called "Fat-Higgs" models [isj^ in which 
some TeV-scale states are composites of the underlying strong-coupling dynamics. This 

^ Other works which consider aspects of dark matter phenomenology in the context of Hidden Valley 
or Higgs Portal models are contained in Ref. [HI [HI, while earHer related studies are contained in Ref. [TU] , 
^ Other models in a similar class to the Fat Higgs model are discussed in Ref. [I9j . 
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UV completion is consistent with both colhder constraints and aesthetic requirements 
such as gauge couphng unification. This completion is discussed in detail in Section [HI 

Furthermore, the existence of partially hidden (secluded) sectors is common in models 
that attempt to embed the SM within a larger structure. Well studied examples include 
higher-rank GUT models, such as those based upon Eq [20], and supersymmetry breaking 
models, in particular the messenger sectors of gauge-mediated supersymmetry breaking 
models [5]. More recently, it has been argued that secluded or hidden sectors in the form of 
Randall-Sundrum-like warped "throats" are a ubiquitous feature of the landscape of 
string compactifications [22j, thus implying that there is not an insignificant probability 
that a hidden or secluded throat with a mass scale close to the weak scale exists. In 
fact, as argued by Patt and Wilczek [3], the scales in sectors interacting by Higgs portal 
interactions are commonly tied together. 

Naturally, if our dark matter candidate is to be the dominant component of the cos- 
mological dark matter, we must ensure that the usual neutralino dark matter candidate 
of the MSSM leads either to a subdominant relic density or is unstable. In the case 
in which i?-parity is conserved and a neutralino is the lightest supersymmetric particle 
(LSP), the thermally generated abundance of such a state is in many models well below 
the measured dark matter density. In particular, wino-like or higgsino-like LSPs annihi- 
late very efficiently, leading to subdominant abundances [23] • Coannihilations with other 
supersymmetric states can also deplete the neutralino abundance in many models |24j . 
Alternatively, instead of being a neutralino, the LSP could be a different supersymmetric 
state, such as a gravitino. Within the context of gauge-mediated supersymmetry break- 
ing, for example, the LSP is typically a light gravitino which constitutes only a very 
small fraction of the cosmological dark matter abundance. On the other hand, if there 
exist i?-parity violating interactions, then the LSP will be unstable thus evading this 
issue entirelylfl 

Turning to the structure of our paper, in Section [2] we introduce our models and 
explain how they are a modified form of the so-called Minimal Non-minimal Supersym- 
metric Standard Model (MNSSM), while in Section [3] we give a brief introduction to 
the physics of the Sommerfeld enhancement that plays an important role in our calcula- 
tions. In Section m we summarize the calculation of the relevant dark matter annihilation 
cross section including the Sommerfeld enhancement and present our results for the relic 
density. In Section [5] we briefiy discuss the direct and indirect detection of our dark 
matter candidate, leaving a more detailed study for a companion paper [T7]. Sectional 
in which we demonstrate that our models may be given a UV completion in so-called 
"Fat-Higgs" models where the states are composites of underlying strongly coupled dy- 
namics, is somewhat outside the main development of our paper and may be skipped by 
readers only interested in dark matter phenomenology. Finally, our conclusions are given 
in Section [3 

■^A late-decaying LSP may even be beneficial in that it can correct the BBN prediction for the ^Li 
to '^Li ratio [25]. 
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2 The Supersymmetric Higgs Portal Model 



The relevant terms of the model that we wish to study are specified by the superpotential 

W = Wmssm{i^ = 0) + XNH^H, + -NS' + -^S' + (1) 

where N and 5* are SM singlets and gets an electroweak-sized scalar vacuum expecta- 
tion value (vev). The term Wms5m(/^ = 0) refers to the MSSM superpotential without 
the "/i" term, while the elhpsis denote terms, such as possible tadpoles, that will not be 
important. 5* has an exact non-R Z2 symmetry which will be unbroken in the vacuum 
and which leads to a stable relic, s, the fermionic component of the 5* superfield with 
mass rris. Note that N does not have a mass term before electroweak symmetry is broken. 
In fact the masslessness of before EWSB is not crucial; all that is required is that 
the mass of is small compared to S as we explain in detail below. We will assume 
that the standard neutralino supersymmetric dark matter candidate is irrelevant, either 
because Rp is broken, or because its relic density is subdominant to that of s. This model 
Eq-dID is a simple variation of the model outlined in Refs. [27] and [28], referred to as the 
Minimal Non-minimal Supersymmetric Standard Model (MNSSM) in Ref. j^, where 
the superpotential has the form Wmnssm = WMssAiif^ = 0) + XNHuHa + ^2-^, and ^2 
is a mass dimension two "tadpole" -term parameter that is in general possible. 

The annihilation cross sections determining the number density of our dark matter 
particle will depend on the couplings, A and A'. As we will argue in later sections, the 
most interesting dark matter phenomenology occurs when the coupling A' is large. Fur- 
thermore, for large A there are important additional contributions to the higgs quartic- 
self-couplings and the upper bound on the lightest higgs mass is considerably raised. 

For sufficiently large A, A' couplings the theory Eq.([T]) hits a Landau pole below the 
Planck-scale, and so must be considered a low-energy effective theory with a cutoff A. 
We will argue in Section [6] that the above effective theory can result from a limit of 
the Fat Higgs model [18] in which S" is a composite meson field of new supersymmetry- 
preserving strong-interaction dynamics, giving our effective theory a possible and plau- 
sible UV completion without tadpole problems, and also with a natural reason to expect 
large couplings A and A'. We emphasize that the Fat Higgs model is only one of many 
possible examples of a UV completion with large couplings A and A'. To analyse the 
dark matter phenomenology it is sufficient to focus on the effective superpotential in 
Eq.([T]) rather than that of any particular UV completion. Although the precise form of 
the superpotential, Eq.([T]), arises as a low-energy limit of the Fat Higgs model, other 
terms, such as a bare //-term, a mass term for A^, or A^^ self-interactions can be added to 
Eq.([T]) without qualitatively changing our results if the final mass of A^ is parametrically 
smaller than by a factor of (9(10). We assume this in the remainder of our analysis. 

From the superpotential, Eq.([T]), the Lagrangian terms determining the important 

^More precisely, the spectrum of neutral Z2-even Higgs scalars arising from i7„, Hd, N after mixing 
must be such that a state with substantial interaction with the Z2-odd states s, s has mass which is 
parametrically small compared to m^. 
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interactions and masses in the model are 



-^fermion -^scalar; (2) 



where, 



~ ~ ~ ~ A' Tfi~ 

Acrm = —^nhuhd — Xnhuhd — Xnhuhd — —nss — X'hss ^ss + h.c. + 



Cscai = \X'ns + mss\^ + 



X' ^ 
XKhd + —s^ 



+ \Xnhd + XttLtR\^ + \Xnhy? 



+ soft susy breaking terms + 



where n (n), s (s), (hu) and hd (hd) are the scalar (fermionic) components of the 
superfields A^, S, and Hd respectively. 

To simplify the analysis, we make the reasonable assumption that the scalar tri-linear 
A-terms and bilinear B-terms are small and consequently we neglect their effects in cross 
sections. In particular we are setting the tri-linear A-terms, Ax = Ax/ = 0. We also 
neglect the D-term interactions as these give irrelevant 4-point Higgs interactions. 

To assess the viability of our dark matter candidate, we need to calculate its thermal 
relic abundance. An important point to note is that the freeze-out temperature of our 
dark matter particle is higher than the electroweak phase transition temperature, Tc, 
for the range of dark matter masses ^ 3 TeV we consider. (In our companion paper 
[T7] we will explore the region of dark matter masses below 3 TeV.) Consequently, in the 
relic density calculation, electroweak symmetry is still a good and no Higgs scalars will 
have vevs. Moreover, above Tc, the fermionic states n, hu and hd are massless, as are all 
quarks and gauge bosons. The only massive fermionic state of interest is s with mass mg. 
In the scalar sector, the thermally-corrected masses of the scalar n states and MSSM 
Higgs bosons are taken to be negligible compared to m^, which is a good approximation 
for the parameter range we are interested in. 

Taking rris ^ 3 TeV does lead to one slight complication in our analysis in that the 
scalar state, s, has a very similar Boltzman factor compared to s near the freeze-out 
temperature, Tf. This is due to the fact that the mass splitting between s and s is small 

rris - m-s = {ml + ml^^^f^ - ~ ml^^^/m-, < Tf ~ m,-/25, (3) 

where msusy is the supersymmetry breaking scale, which is parametrically smaller than 
rris. This means that the scalar s and fermion s states will freeze-out at roughly the same 
temperature and we have to consider the annihilation rates of the scalar states as well 
as the fermionic state^. 

In addition to the purely scalar interactions which follow from Eq.Q the fermionic 
interactions which are of importance in determining the relic abundance of our dark 



^We remark in passing that our qualitative conclusions regarding the dark matter freeze out density 
would not be changed if a scalar component of S were the lightest Z2-odd state, say due to CP- violation. 
The Sommerfeld effect acts equally for both scalar and fermionic annihilating particles as explained in 
Section 3. 
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matter state are 

- ■^^{X'4>nSMCsM +iX'anSljC'j5SM) 

- KKoPLnM - {KT^m^PrKd - {KYKbPlUm - K^m^^rKd^ 

where Pl,r = (1 ±75)/2 and we have rewritten the fermionic states in terms of Majorana 
and Dirac spinors indicated by the subscripts M and D respectively. The scalar states 
have been written in terms of their CP-odd and CP-even components, denoted generically 
as Ai = + ^ctj), and C is the charge conjugation matrix. The subscripts u and d on 

the scalars refer to the Higgs "up" and "down" states. 



3 The Sommerfeld Enhancement 

For dark matter particles moving at small relative velocities, the exchange of scalar states 
leads to an enhancement by factors depending on the inverse velocity, 1/v. This Sommer- 
feld enhancement corresponds to the summation of a series of ladder diagrams where the 
scalar state is repeatedly exchanged (see Fig. [1]). This enhancement is only significant if 
there exists an S-wave annihilation amplitude, otherwise the angular momentum barrier 
will suppress the effect |§ 




Figure 1: Generic Sommerfeld diagram. The "blob" vertex represents all possible 5- wave 
annihilations of the incoming states including s-channel, t-channel and annihilation via 
contact interactions. 

The calculation of the Sommerfeld enhancement can be formulated in terms of a 
non-relativistic quantum two-body problem with a potential acting between the incoming 
particles. This is equivalent to the distorted Born-wave approximation common in nuclear 

^If vector states are exchanged, there can either be an enhancement or suppression depending on the 
relative charges of the annihilating particles. 
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physics. To a good approximation this leads to a dressing of the S'-wave part of the tree- 
level cross sections with a multiplicative factor, 



^ = ^4?e- (4) 

The full calculation of R can be involved and in many cases, including that of a 
Yukawa potential, cannot be solved analytically. In our model the only particles which 
can act as the "rungs on the ladder" in the Sommerfeld diagram shown in Fig. [1] are the 
scalar n states, the n fermions not contributing to the enhancement. The non-relativistic 
potential which is relevant for all the diagrams we will consider is found to be 

V = -^e-™"^ (5) 

where is the mass of the particle acting as the "rungs on the ladder" . The Schrodinger 
equation for the two dark matter particle state, ip, with this potential reads 



mz dr 



" + \/.^ = K^^ (6) 



where K = Mv"^ is the kinetic energy of the two dark matter particles in the center-of- 
mass frame, where each dark matter particle has velocity v. Using the outgoing boundary 
conditions, '?/''(oo)/V^(oo) = inisV, R is given a.s R = \ip{0)/ip{oo)f. In the simple case 
we are considering, we can derive an analytic form for R. In the limit where the ratio 
e = nin/ms = 0, R takes the form jl6j . 

R= , ^ , (7) 

where y = /8v = /Avr and Vr = 2v is the relative velocity between the two dark 
matter particles. Taking the small Vr limit we have 

A'" 

R-— (8) 
and we see that this effect will be largest for small Vr- 



4 Calculation of the Relic Density 

We are now in a position to calculate the relic density of our dark matter candidate. 
As mentioned in the previous section, in this paper we will restrict our analysis to dark 
matter particles with masses mg > 3 TeV. Not only is this range of masses physically 
interesting, it also simplifies the analysis considerably due to the fact that freeze-out oc- 
curs at a temperature above the electroweak phase transition, thus leading to a situation 
in which no scalars have vevs. Consequently, the number of possible vertices contributing 
to the annihilation cross sections is reduced and the calculation of the relic abundance 
greatly simplified. 
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There are three important types of diagram which determine the rehc abundance of 
our dark matter particle. The first type (type I) involves the annihilation of two scalar s 
states. We can have two scalar s states annihilating into Higgs, Higgsinos, scalar n states 
or fermionic n states as depicted in Fig. [2l For computational ease we take all states to 
be massless apart from s and s which have masses and mg + ml /rris respectively. 




Figure 2: Type I annihilation diagrams for the scalar s states. 



For all scalar annihilation diagrams we receive an enhancement from the Sommerfeld 
effect where the CP-even scalar, acts as the "rungs on the ladder" between the 
annihilating scalar s states as depicted in Fig. [3] for the case of the annihilation of two 
s states. The "blob" vertex represents all possible ways of annihilating the s states 
including s-channel, t-channel and annihilation via the 4-point vertex. 



X 



s) 



Figure 3: Sommerfeld diagram for scalar annihilations. 



The resulting self annihilation cross sections for the CP-even and CP-odd components 
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of s are 



a{(f)s4>s [huhd]) = (j{asas [Khd]) 
a{(j)s(t)s [Khd]) = cr{asas -> [Khd]) 



a{as(l)s {Khd)) 



0- 



nn] 



(A^A)^ y 

327rw^m| 1 - e-s^' 

(A^A)^ y 
327rt>^m? 1 — e~2/' 



(A-A)^ y 



a{(j)s(ps [nn]) = aia^a^ [nn]) 



y 



0" OoO, — i> nn 



1/ 



y 



IQTTVrmj 1 — e y 

where the factor y/{l — e^^) accounts for the Sommerfeld enhancement and 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 



[Khd] = mi.Kh-dAKh-d)\ 

[huhd] = (t)u(t)d,auad,h'![h'^,{K[h';^)* , 

{huhd) = (t)uad,au(t)d,h1;^h~,{h'^h~y, 

[nn] = 0„0„, a„a„, 

represent all possible final states in each case. 

The second type of diagram (type II) we need to include is the annihilation of a 
scalar s state with a fermionic s state. The relevant diagrams are shown in Fig. HI Each 
process in Fig. H] can also be enhanced by the Sommerfeld effect via the diagram shown 
in Fig. [5], where the "blob" vertex represents both s-channel and t-channel processes. 
The enhancement factor is exactly the same in this case as it was for the annihilation of 
scalar s states. 



''h h 

/ iijui"'a 



n 



n 



Is 



hdi hu 



n 



Figure 4: Type II: Annihilation of s with s. 
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Figure 5: Sommerfeld diagrams for annihilation of s and s. 



The cross sections for these processes are found to be 



a{(f)sS (hihj)) = cr{asS {hihj)) 
a{(f)sS — ^ {nh)) = a{asS — > {nh)) = 



(A'A)2 



y 



y 



327rt>j,m| 1 — e^y 



where 



{nh) = anh, ^^n, 



(15) 
(16) 



represent all possible final states for each process. 

The third type of process (type III) is the annihilation of the s states. In the elec- 
troweak symmetric limit, there are no vevs for the n or Higgs states, neither is there a 
tri-linear scalar A term, Axnhuhd as we have approximated this term to be zero. This 
means that the only Higgs-like final states from s annihilation will be products of neu- 
tralinos and charginos. We can also have t-channel exchange of a s, which produces a 
CP-odd final state pair, an0„. All possible diagrams with non-zero S'-wave amplitudes 
are shown in Fig. [61 




Figure 6: Type III: Annihilation of two s states. 

The ss annihilation processes can be enhanced by the Sommerfeld factor if the initial 
s pair are in an S'-wave state. The corresponding Sommerfeld diagram for ss annihilation 
is shown in Fig. [71 where the "blob" vertex represents both s-channel and t-channel 
processes. 
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Figure 7: Sommerfeld diagrams for the annihilation of two s states. 
The resulting cross sections for the s annihilations are found to be 

a{rs^[h.M) = ^ -y ^ (17) 

a{s~s^hh) = ^ (18) 

^ ' 647rf^m| I - e-y ^ ' 

where [/im/z.^] was defined earlier in this section. 

We now have all of the cross sections needed to determine the relic density. As we 
have two states freezing out almost simultaneously (our dark matter state s and its scalar 
partner s) we must be careful to include the effects of the heavier state in the calculation 
of the relic abundance of the dark matter species. We follow Refs. [Ml [29] in calculating 
the final relic abundance of our dark matter candidate. 

If we relabel our two states, s and s, as si and S2 respectively, the type of reaction 
that will determine the freeze-out of our two particles is 



CTij = Cr[SiSj 



XX'), (19) 



where X and X' will be some combination of Higgses, higgsinos, fermionic n states and 
scalar n states, which will decay to lighter MSSM degrees of freedom. Taking into account 
all possible diagrams, the three cross sections we are concerned with have the following 
forms 

a{s,s,-.XX) = 64,,^^. (20) 



a{s,s,^XX) = (21) 



rr(. . ^ XX'\ [3(A-A)^ + 7{Xr] y .... 

where we have averaged over the components of the initial scalar states where appropri- 
ate. 

We assume that any S2 states remaining after freeze-out will eventually decay down 
to siXX'. This means that the total number density of our dark matter particle will be 
equal to the sum of the Si and S2 number densities at freeze-out. 



10 



In order to calculate the relic density we define the following useful quantities [23] 



where 
and 



(7.(l + A0'/'exp[-xA,] 
fl'cfr 



Ai = {rrii - rriij/mi, 



9es = ^9i{'^ + Ai)^/^ exp[-xAi 

i=l 



(23) 

(24) 
(25) 



where Qi is the number of degrees of freedom of Sj, is the mass of Sj and x = rris/T. 
Of course in our case we only have two different species of particle and so only A2 is 



non-zero. In fact as Si = s and S2 = s, we have A2 



/rris. Each of 



our Si states have gi = 2 degrees of freedom. Following Ref. |24j, we find the freeze-out 
temperature, T/, by iteratively solving the equation 



X/ = In 



0.038gesMpims {cTeffVr) 

1/2 1/2 
9* X/ 



(26) 



where xj = rris/Tf and 



Ceff 



2J 



a,j^{l + A,;)3/2(i + A,)'/' exp(-x(A, + A,)). 

fl'cfT 



(27) 



The final relic density is given by 



1.07 X lO^Xf 

GeV)j' 



where J 



X Oeff 



dx 



{2i 



and Qi, is the total number of relativistic degrees of freedom at Tj. In our calculation 
of the relic density we will take (7* = 248, which includes all MSSM degrees of freedom 
plus the four associated with the extra superfield A^. In order for this to be correct the 
masses of all these states must be below Tf ~ ms/^b, which will be true when we take 
rus > 3 TeV and msusy = 100 GeV as an example parameter set. 

It is instructive to compare the two cases of when we correctly include the Sommerfeld 
factor in cross sections and when this contribution is absent. The comparison is most 
clear when we plot the relic density, fig/i^, against as shown in Fig. [HI In Fig. [SI 
the red dashed lines correspond to the case where the Sommerfeld factor (given by 
R = y/ (1 — e~^)) is not included in the cross sections, while the blue solid lines correspond 
to the case where it is. The three lines for each case (with and without the Sommerfeld 
effect), starting from the furthest left, correspond to A = A' =2, 2.5, 3 respectively. 
The two lines parallel with the axis correspond to the WMAP allowed range for the 
dark matter relic abundance, inferred from the combination of QmH"^ 
VL^h'^ = 0.02229 ± 0.00073 [30j. 



0.12771™ and 
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Figure 8: flh'^ as function of at fixed A and A'. The red daslied lines correspond to 
the case where the Sommerfeld correction is not included where as the blue solid lines 
correspond to the case when it is included. The furthest most left line for each colour 
corresponds to when A = A' = 2, the middle lines are when A = A' = 2.5 and the lines 
furthest right are when A = A' = 3. All plots are produced using mgusy = 100 GeV. 

For each line (of fixed coupling), the relic density increases as we increase the mass, 
rris, as we would expect. Comparing sets of contours with the same couplings (A = 
A'), we see the dramatic effect of the Sommerfeld enhancement. When the Sommerfeld 
enhancement is included, the annihilation cross sections are increased, thus depleting the 
number density of the dark matter particles which survive after freeze-out. The bottom 
line is that the Sommerfeld enhancement allows for very heavy dark matter particles 
to provide the required dark matter relic abundance. From Fig. [8] we can see that the 
maximum mass consistent with the WMAP allowed range when we have A = A' =3 is 
close to 25 TeV. 

The results of a full numerical scan (including the Sommerfeld enhancement) over the 
three parameters A, A' and mg is shown in Fig. [9l Here, we plot contours corresponding 
to the allowed range of Qsh"^ in the A — A' plane. Each pair of contours correspond to a 
different value of the mass, mg, between 3 and 23 TeV. The left (right) contour of each 
pair corresponds to the higher (lower) end of the allowed range in Qgh"^- 

Although we show contours only for discrete choices of ms, the remaining regions of 
the A — A' plane are filled for intermediate values of the dark matter massjj The effect 
of the Sommerfeld enhancement is to pull the pairs of contours downward towards the 
bottom left corner of the A — A' plane. This allows us to have the correct relic density 
for a given dark matter mass for smaller values of the couplings. 

On examination of the parameter region shown in Fig. [9|, it is interesting and perhaps 

^There will be an upper limit on how large the couplings can be, which is determined by insisting we 
have perturbativity up to our cut off scale. As shown in Section [SI large couplings of the size considered 
here are shown to be natural in a consistent UV completion. 
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!j= 11 TeV 
ij= 15 TeV 
i~= 19 TeV 
23 TeV 

S 



1.5 2.0 2.5 3.0 



Figure 9: Plots of pairs of contours for the allowed range of Qh"^ in the A — A' parameter 
plane for different values of the mass, m^. We have contours corresponding to masses 
3, 7, 11, 15, 19, 23 TeV. The contours are produced using msusy = 100 GeV. 

important to note that we are able to generate the correct dark matter relic density 
(for a given mass) for relatively small A couplings, provided we have a large enough A' 
coupling. The main reason for this is that only the A' coupling appears in the Sommerfeld 
enhancement and as we can see from Fig. [HI it is the Sommerfeld enhancement that 
allows us to have large masses for the dark matter particles. With this in mind, it is also 
interesting to note that as we have annihilation diagrams which depend on A' only, we 
can take A to be small (around 0.1 for example) and still have viable dark matter with 
masses up to around 23 TeV. 



5 Direct and Indirect Detection 



Although we intend to discuss the prospects for the direct and indirect detection of 
heavy Higgs Portal dark matter in some detail in a companion paper to follow [17] we 
will here briefly touch upon this subject. We find that the direct detection phenomenology 
is fairly conventional and, although present experiments do not yet lead to restrictive 
limits, a sizeable fraction of the expected parameter space will be covered by proposed 
next generation detectors. In contrast, the indirect signals are greatly modified by the 
potentially very large Sommerfeld enhancements. 
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5.1 Direct Detection 



Experimental programs designed to observe the elastic scattering of WIMPs with nuclei 
are collectively known as direct detection. The dark matter particles in our model, s, 
interact with quarks in nuclei through the effective scalar interaction given by 

C= Cu~s~sUU+ CdssDD, (29) 

U=u,c,t D=d,s,b 



where 



Cu = Y.—^--'^ = Y.^^^ (30) 



and the mixing matrix Vij specifies the admixture of n, and h^^ states in the neutral 
scalar mass eigenstates, hi, with lightest neutral Higgs state being denoted hi. Unlike in 
the case of many other dark matter candidates, there is no contribution from Z exchange 
in this model. Note that the kinematics of the interaction (even if we consider scattering 
of individual nucleons or even quarks with the dark matter) are such that we are outside 
of the range for which the Sommerfeld enhancement is important. 

Following Refs. [SUES], we estimate that this interaction leads to an elastic scattering 
cross secton per nucleon of 



7 . ^ViA" fX'V /120GeVA' , , 



For the range of masses we are interested in here, this cross section is below the 
current constraints from experiments such as XENON [52] and CDMS but is likely 
to be reached in the next few years. For less optimal values of A', uih^ or Vij, however, 
the prospects for direct detection could be considerably more difficult. 



5.2 Indirect Detection 

In addition to direct searches for dark matter, astronomers are also searching for the 
products of dark matter annihilations, including gamma-rays, neutrinos, positrons and 
antiprotons [34j. These efforts are known as indirect detection. 

The dark matter annihilation rate, and thus indirect detection rates, can be enor- 
mously enhanced due to the Sommerfeld effect. Depending on the astrophysical envi- 
ronment being considered, annihilation rates can be enhanced by factors of 10^ to 10^ 
or even greater due to the slow relative velocities of dark matter particles. In fact, the 
velocity dependence of the enhancement factor can potentially favour such astrophysical 
objects as dwarf satellite galaxies of the Milky Way (due to the extremely low velocity 
dispersion) as sites for indirect detection, rather than the central regions of the Milky 
Way itself. A full calculation of the expected flux depends upon a detailed knowledge 
both of the resonance structure of the Sommerfeld enhancement in the non-coulombic 
and low Vr regime and of the sizes of vacuum expectation values and interaction terms 
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in the scalar [S, A^)-Higgs sector. A preliminary estimate shows that current indirect 
detection experiments do not impose a useful limit on heavy Higgs portal dark matter, 
but that there is a potential for significant signals in future observations [17j . 



6 UV Completion as a Fat Higgs Model 

It should be noted that the sizes of the couplings we have taken in the analysis of 
section H] are the values for typical momentum transfers at freeze-out (~ m^/^jo, with 
Pfo ~ 0.2). The cut-off of our effective theory will be related to the energy scale at which 
our couplings become non-perturbative, which for definiteness we take to be where the 
two loop terms in the renormalisation group equations for A and A' become of order the 
one loop terms. For example, if we take = 3TeV, A = 0.8 and A' = 1.6 the cut-off 
is ~ 4000 TeV. A more extreme example is where we take rris = 23TeV, A = 1.0 and 
A' = 3.2 with cut-off ~ 70 TeV. The consequence of having a cut-off below the GUT 
scale is that we are motivated to think about how this model can be UV completed. 
We emphasise that for the analysis of the dark matter properties and thermal freeze out 
the effective low energy lagrangian, Eq.([T]), is appropriate as there is a large separation 
between mg/3/o and the cut-off even for the most extreme case we consider, = 23 TeV. 

One possible way to UV complete our model and justify the choice of large couplings 
A, A' is to have some strongly interacting physics which dynamically generates the su- 
perpotential S mass. It is noteworthy that the "Fat Higgs model" of Ref. [18] provides 
exactly such a UV completion. With this in mind, we will now describe how our effective 
theory can arise in a certain limit of the Fat Higgs models. 

The Fat Higgs model is an = 1 supersymmetric SU{2) gauge theory with six 
doublets with the quantum numbers shown in Table [H 

The tree-level superpotential is given a^ WpHtot = Wi + W2 + W3 where 



The P and Q mixing terms are there to marry off unwanted "spectator" states such 
that the low energy effective theory is as minimal as possible. It is also possible to apply 
a Z3 which protects us from tadpole terms involving either of the singlet fields. Si and 
5*2 [18]. This will commute with the existing symmetries. 

The gauge symmetry SU{2)h becomes strongly coupled at some scale, Kh- Below A//, 
the appropriate degrees of freedom are mesons which are composite objects consisting 
of two "T" doublets in the form Mij = TiTj, with (i, j=1...6). There is a dynamically 
generated superpotential of the form PfM/A^ as well as the tree level superpotential 

®The terms with coefficients j/3 and 1/4 were not included in Ref. [18 . These terms are not forbidden 
by any symmetries so we include them for completeness. 



W2 



yiSiTiT2 + y2S2TsT4 + y^SiT^T^ + y4,S2TiT2 
-mT^TQ 



(32) 
(33) 




(34) 
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Table 1: The field content under an S\J{2)l x SU{2)h gauge and ^f/(2)^j x SU{2)g x U{\)n 
global symmetries. There is also an accidental Z2 symmetry with fields transforming as 
shown. The f/(l)y subgroup of SU{i)R is gauged. 



which follows from Eq.( l34l) . As P,Q,Si and 5*2 are not charged under SU{2)h, they 
remain fundamental below Ah- The canonically normalised effective superpotential reads 

Wdyn = A (PfM - v^oM^e) + m^SiMi2 + m2^2M34 + m^^SiM^i^ + miS2M^2 
+ ms (M15P11 + M16P12 + M25P21 + M26P22) 

+ me (MssQn + M36Q12 + M45Q21 + M36Q22) , (35) 

where, using Naive Dimensional Analysis (NDA) [35j, we have 

2 , . 

~ (4^0^' ^''^ 

mi ~ (37) 
X{Ah) ~ 47r. (38) 

We now make the assumption that (777.5, me) 3> (mi, m2, ma, m4), by a factor of 10 
or so, and integrate out everything with a mass proportional to ms or rriQ. This leaves 
us with a superpotential of the form 

W;,y^ = AM56 (M14M23 - M24M13 - t;^ + M12M34) 

+ miSiMi2 + m2S2Ms4 + m3S'iM34 + m4S'2Mi2. (39) 

Assuming that mi ~ m2 ~ 777-3 ~ ^4 ~ ""^'j the fermionic components of the super- 
fields, 5*1, 5*2, M12 and M34, mix and, provided mim2 7^ m3m4, the lightest eigenvalue of 
this mass matrix will generically have a mass of order m'. 
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If we now do this diagonalization and integrate out all but the lightest eigenvalue, 
call it S, of the 5*1, 5*2, M12, M34 mass matrix, we are left with the superpotential 



W = XN {H^H, - vl) + ^NS' + (40) 



where we have changed notation according to the identifications 

H+ \ ( Ml, \ ( m\ ( M, 



_ 13 ^ = 



N = M56. (41) 



The parameter A' = XUijUki, where U-ijUki are components of the unitary matrix that 
diagonalizes the Si, 5*2, M12, M34 fermion mass matrix. The indices on the Us are there 
for show, the basis is irrelevant as we do not really care about the specific mixing between 
states. 

The final assumption we make is that is parametrically larger than the electroweak 
scale and soft supersymmetry breaking masses. The superpotential in Eq.l HTl) is of the 
form we need with an additional linear term for the superfield N . This term is harmless 
with respect to the dark matter dynamics but we include it for completeness. 

Ignoring the S field for now, the remaining superpotential is that of the Fat Higgs 
model and the analysis of the electroweak vacuum structure proceeds as outlined in 
Ref. [IB]. It is worth comparing the superpotential in Eq. (l4ip with that of the MNSSM 
[271 ES]- In particular, the N linear term in Eq. fl4ip is analogous to the tadpole terms 
appearing in the superpotential of Eq.(3.1) of Ref. [27]. In fact, the superpotential in 
Eq. pT]) (apart from the S terms) is that of the MNSSM. Consequently we can use the 
rather more detailed analysis of Refs. [271 128] for the Higgs sector. 



The S terms in Eq. (l4Ti) do not spoil the electroweak structure of the MNSSM. We 
can see this by integrating out 5* using the equations of motion 

dW 

—- = }^NS + m,S = Q, S = 0. (42) 

ob 

Substituting the solution back into Eq. (HTl) we have the effective superpotential 

W,s = XN {HuHa - vl) , (43) 
which is exactly the same as the superpotential for the MNSSM and the Fat Higgs model. 



7 Conclusions 

In this article, we have discussed models in which a very heavy (3-30 TeV) dark matter 
candidate is present. In particular, we have focused on models motivated by Higgs Portal 
and Hidden Valley models, in which the dark matter (and the rest of the partially hidden 
sector) interacts with the Standard Model and its superpartners only through Higgs 
interactions. 
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Dark matter annihilations in this scenario are considerably enhanced by non- 
perturbative contributions known as the Sommerfeld effect. Through this enhancement, 
dark matter particles with masses well above the electroweak scale can be produced ther- 
mally in the early universe with an abundance consistent with the measured density of 
dark matter. The dark matter particle in this scenario, although well beyond the reach 
of the Large Hadron Collider, is still potentially detectable by direct and indirect dark 
matter experiments. Although we leave the details of this to future work [L7\, we point 
out that Sommerfeld corrections can dramatically enhance the dark matter annihilation 
rate in low velocity dispersion environments, such as dwarf spheriodal galaxies, thus 
considerably improving the prospects for indirect dark matter searches. 

The particular model we study, which adds two extra SM singlet states to the MSSM 
spectrum, is independently motivated by the desire to raise the upper bound on the 
lightest higgs mass, thus lessening the LEP fine-tuning constraints. We also showed that 
our model may be given a UV completion in the form of the previously considered "Fat 
Higgs" models, where the singlet states are composites arising from strong hidden-sector 
dynamics. 

Thus, Higgs Portal Dark Matter provides an example of an attractive and motivated 
alternative to conventional MSSM neutralino dark matter which is less fine-tuned and 
may be tested by current and future indirect detection experiments. 
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